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IONOSPHERIC  EQUIVALENT  SLAB  THICKNESS 
AND  ITS  MODELING  APPLICATIONS 

1.  INTRODUCTION 

There  are  a  variety  of  ionospheric  models  currently  in  use  or  under  development, 
each  with  its  own  approach  to  the  problem  of  describing  the  global  ionosphere.  Theoret  ical 
models  have  grown  greatly  in  sophistication  over  the  course  of  this  decade,  both  because 
of  the  vast  increase  in  computing  power,  and  because  of  the  growth  through  data  analyses 
in  our  qualitative  understanding  of  the  processes  that  occur  in  the  ionosphere.  One  of 
the  more  well-known  of  these  theoretical  models  was  developed  at  Utah  State  University 
(e.g.  Sojka  and  Schunk,  1985),  where  a  time-dependent,  three-dimensional,  multi-ion 
numerical  solution  is  achieved  for  the  global  F  region.  This  model  achieves  good  qualitative 
descriptions  of  the  ionosphere  (for  example,  in  the  equatorial  anomaly  region,  the  mid¬ 
latitude  trough,  and  the  winter  anomaly).  However,  the  quantitative  agreement  is  not 
always  good.  While  this  model  is  being  further  validated  and  refined  against  now  or  existing 
data  sets,  and  while  these  numerical  experiments  are  important  to  our  understanding 
of  the  global  ionospheric  morphology,  the  general  scientific  community  continues  to 
work  with  empirical  or  semi-empirical  ionospheric  models.  'Phis  is  not  only  due  to  their 
agreement  with  known  data,  but  also  for  ease  of  numerical  computation.  And,  in  many 
cases,  specification  of  the  mean  (or  quiet  condition)  ionosphere  is  quite  adequate. 

The  basic  difference  in  emphasis  between  theoretical  and  semi-emp  .  ical  models  is 
that  the  former  seek  a  self-consistent  global  solution,  while  the  latter  define  the  electron 
density  profile  at  any  given  location  and  under  any  conditions.  One  such  semi-empirical 
model  is  the  International  Reference  Ionosphere,  or  IRI  (Re.vcr,  1981).  This  combines 
empirical  global  maps  of  the  profile  peak  parameters  with  analytic  expressions  for  the 
shapes  of  each  region  of  the  ionosphere  for  the  desired  location,  season  and  solar  activity 
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conditions.  The  Bent  (Llewellyn  and  Bent,  1973)  model  underwent  a  similar  development, 
but  utilized  a  different  set  of  global  maps  and  relied  on  simpler  expressions  for  the  profile 
shapes.  One  model  under  continued  active  development  is  the  US  Air  Weather  Service 
Ionospheric  Conductivity  and  Electron  Density,  or  ICED,  model  (Tascione  et  al.,  1988). 
This  model  seeks  to  incorporate  various  processes  or  features  of  specific  regions  of  the  globe 
(such  as  the  SLIM  model,  Anderson  et  al.  (1987),  at  low  latitudes)  into  a  smoothly  varying 
global  model,  and  it  is  also  designed  to  respond  to  real-time  ground-based  and  satellite 
observations.  This  latter  point,  making  use  of  all  available  real-time  observations,  is  a  very 
positive  step  in  the  field  of  ionospheric  modeling. 

For  the  reasons  stated  above,  in  the  majority  of  ionospheric  research  today, 
specification  of  the  electron  density  profile  is  the  key,  as  this  will  drive  the  basic  ionospheric 
model.  The  peak  parameters  (such  as  foF2  and  h’F2)  are  generally  obtained  from 
empirical  or  semi-empirical  global  maps.  These  are  simple  numerical  representations  of 
observations  and  can,  in  principle,  be  constructed  to  any  desired  accuracy  over  a  particular 
sounding  network,  according  to  the  particular  problem  being  addressed.  The  profile 
shape,  however,  is  a  complicated  function  of  the  neutral  and  plasma  temperatures,  ionic 
compositions,  neutral  winds,  and  therefore  of  the  season,  local  time,  location,  level  of  solar 
and  geomagnetic  activity.  Qualitatively,  the  effects  of  varying  each  of  these  parameters 
on  the  profile  shape  is  reasonably  well  understood,  but  the  problem  of  obtaining  a  reliable 
physical  (or  numerical)  specification  of  the  variations  has  not  been  achieved  to  date.  There 
is  a  continued  process  of  model  validation  and  refinement. 

The  equivalent  slab  thickness  is  defined  as  the  breadth  (in  km)  of  an  ionosphere  of 
uniform  electron  density  equal  to  the  peak  density,  Nmax.  Functionally,  r  =  TEC/Nmax, 

OC 

where  I  EC  =  /  Nc(h)  dh,  and  it  is  a  convenient  one-parameter  summary  of  the  election 

o 

density  profile.  Indeed,  it  can  be  related  to  a  variety  of  quantities  of  interest  that  effect 
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the  overall  profile  shape.  Apart  from  the  intuitive  relation  to  the  F2  layer  semi-thickness, 
various  authors  have  proposed  relations  between  r  and  both  neutral  temperature  and 
0+ / H+  transition  height  (Titheridge,  1973),  the  scale  height  of  atomic  oxygen  and  drift 
motion  (Furman  and  Prasad,  1973),  and  the  mean  gradient  of  electron  temperature 
(Amayenc  et  ah,  1971).  In  addition,  r  is  operationally  a  very  useful  parameter,  as  it  allows 
one  to  make  a  simple  conversion  between  foF2  and  TEC. 

The  equivalent  slab  thickness  is  a  parameter  with  largely  unexplored  applications  in 
the  field  of  ionospheric  modeling.  Firstly,  it  is  a  new  database  that  should  be  used  for  model 
validation.  Secondly,  there  is  a  different  emphasis  in  studying  r.  For  example,  validations 
based  on  foF2  indicate  the  models  that  successfully  describe  the  layer  peaks,  but  yield  no 
information  on  the  heights  or  the  general  profile  shape.  Validations  that  use  observations 
of  TEC  largely  assess  the  topside  profile  as  anchored  to  the  layer  peak.  While  r  may  not  be 
easily  interpreted  in  terms  of  neutral  temperature,  winds,  or  composition  changes  (though 
each  of  these  effects  may  be  anticipated),  the  variations  are  a  good  indication  of  how  the 
broad  structure  of  the  electron  density  profile  changes  as  a  function  of  season,  time  of  day 
and  solar  activity.  Thus,  studying  r  does  constitute  a  different  way  of  utilizing  the  available 
databases  of  TEC  and  Nmax. 

Describing  the  observed  variations  of  r  at  a  reliable  site  is  thus  a  worthwhile  exercise, 
and  this  forms  the  core  of  this  report.  In  subsequent  sections,  the  database  is  introduced, 
the  variations  of  r  are  described  both  qualitatively  and  numerically,  and  modeling 
applications  (validation  and  adjustments)  are  explored. 


2.  THE  DATABASE 


The  bulk  of  the  analysis  contained  in  this  report  comes  from  mid-latitudes  in  the 
American  sector.  Measurements  of  TEC  were  obtained  from  Hamilton,  MA,  by  observing 
the  Faraday  rotation  of  signals  from  the  satellites  ATS-3  and  ATS-5.  The  420  km  sub- 
ionospheric  point  of  this  raypath  lies  at  nearly  the  same  latitude  as  an  ionosonde  located  at 
Wallops  Island,  VA,  but  shifted  by  twenty  minutes  of  local  time. 

Hourly  values  of  r  have  been  calculated  using  the  observed  TEC  (shifted  by  15 
minutes  to  be  the  LT  equivalent  over  Wallops  Island)  and  foF2  values.  The  sample  of  data 
consists  of  sparse  observations  from  1965  and  1967,  and  nearly  complete  observations  from 
196S-19S6.  This  gives  us  access  to  two  solar  cycles  of  TEC,  Nmax  and  r.  Possibly  just 
as  significant  is  the  fact  that  the  first  cycle  was  of  average  strength,  while  the  second  was 
the  second  largest  on  record.  Another  significant  point  is  that  the  TEC  values  have  been 
verified  by  experienced  personnel  at  AFGL,  and  should  be  regarded  as  reliable. 

An  important  point  to  bear  in  mind  is  that  slab  thickness  data  combines  the 
errors  of  two  observational  parameters,  TEC  and  Nmax,  each  of  which  has  some  error  of 
measurement.  Also,  each  quantity  is  capable  of  interpretive  error  in  the  measurement. 
Misidentified  layers  on  an  ionogram  will  lead  to  an  incorrect  foF2  (and  hence  Nmax),  and 
an  incorrect  identification  of  the  n7r  ambiguity  (whether  on  a  strip  chart,  or  by  an  automated 
data  analysis  routine)  will  mean  systematically  incorrect  TEC  values.  In  addition,  there 
are  physical  differences  in  the  measurements  of  foF2  and  TEC.  The  former  is  deduced  from 
vertical  (or  nearly)  incidence  on  the  ionosphere,  while  the  latter  is  deduced  from  a  slant 
path.  The  physical  difference  between  the  exact  sub-ionospheric  point  and  the  location  of 
the  ionosonde  will  mean  that  TEC  and  Nmax  will  not  necessarily  relate  to  the  same  profile 
when  horizontal  gradients  are  large,  such  as  in  the  period  immediately  following  sunrise. 

For  all  the  above  reasons,  we  expect  the  database  of  hourly  r  values  to  be  quite 
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noisy.  Figure  1,  showing  the  overplots  of  daily  values  of  each  parameter  for  a  selected 
month,  demonstrates  this  by  showing  the  variation  of  r  to  be  the  least  smooth.  Because 
of  the  noisiness  of  this  database  we  have  adopted  a  two-pronged  attack  to  the  problem  of 
describing  the  variations  of  r.  Firstly,  the  mean  variations  are  investigated.  A  reliable  and 
representative  mean  is  determined,  dependencies  on  observed  quantities  are  sought,  and  the 
variations  are  described,  both  qualitatively  and  numerically,  through  a  simple  coefficient- 
based  model.  It  is  this  mean  r  behavior  that  has  immediate  modeling  applications,  as 
the  existing  ionospheric  models  are  only  concerned  with  mean  behavior.  Secondly,  some 
indications  of  the  day-to-day  variability  are  sought.  Apart  from  a  simple  statistical 
analysis,  and  a  correlation  analysis  between  TEC,  Nmax  and  r,  we  also  describe  some  of 
the  systematic  trends  seen  in  r  in  the  periods  immediately  following  geomagnetic  storms, 
considered  to  be  the  extreme  cases  of  day-to-day  variability  (see  e.g.  Mendillo,  1978). 

3.  DESCRIBING  THE  VARIATIONS  OF  SLAB  THICKNESS 

This  section  of  the  report  is  divided  into  two  subsections,  concerning  the  mean  and 
day-to-day  variations  of  the  equivalent  slab  thickness  parameter,  t. 

3.1  The  Mean  Variations 

As  described  above,  the  database  of  hourly  slab  thickness  values  is  noisy  (more  so  than 
each  of  TEC  or  Nmax  individually),  and  thus  an  appropriate  mean  needs  to  be  carefully 
defined.  As  is  usual  practice  in  ionospheric  work,  we  are  concerned  with  monthly  mean 
behavior  (the  time  scale  in  the  existing  semi-empirical  models).  How  then  do  we  define  a 
‘mean’  for  each  hourly  value  over  the  course  of  the  month  that  is  a  true  representation  of 
the  data?  The  desireable  criteria  for  such  a  mean  are  that  it  varies  smoothly  over  the  24 
hourly  values  for  each  month,  that  is  is  not  unduly  influenced  by  individual  extreme  (and 
possibly  erroneous)  values,  but  that  it  does  reflect  general  trends  in  the  hourly  values.  The 
arithmetic  mean  was  tested  first,  and  was  seen  to  be  strongly  influenced  by  outlying  values, 
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leading  in  turn  to  a  jagged  daily  variation.  In  an  attempt  to  remove  the  extreme  cases 
of  day-to-day  variation,  an  arithmetic  mean  was  also  evaluated  for  only  those  days  that 
were  geomagnetically  quiet  (Ap  <  10).  However,  this  quiet  mean  was  still  influenced  by 
outliers.  The  median  values  were  more  stable,  but  were  not  influenced  by  extended  tails  of 
the  distribution  of  hourly  values.  The  most  successful  mean  tested  was  the  trimmed  mean. 
This  is  simply  the  arithmetic  mean  of  the  central  50%  of  the  values  in  the  distribution. 
Visually,  this  parameter  showed  the  smoothest  overall  variation,  and  the  most  appropriate 
response  to  different  distributions  of  hourly  values. 

Having  now  defined  the  mean  values,  it  is  a  good  place  to  describe  the  apparent 
variations  of  r  and  to  put  them  into  perspective.  Figure  2  shows  a  three-dimensional 
representation  of  the  trimmed  mean  values  of  r  for  typical  low  and  high  solar  activity  years. 
The  diurnal  variation  (good  examples  of  which  are  shown  in  Figure  3)  is  clearly  different 
from  that  of  TEC  and  Nmax.  Indeed,  TEC  and  Nmax  are  quite  highly  correlated  (sec  the 
next  subsection)  and  r  therefore  highlights  the  differences.  Broadly  speaking,  the  nighttime 
r  values  are  far  less  variable  (in  terms  of  an  annual  variation)  than  the  daytime  values. 

In  particular,  daytime  r  values  are  higher  in  summer  than  winter,  with  summer  daytime 
values  being  higher  than  the  nighttime  values  and  the  reverse  being  the  case  in  winter.  The 
diurnal  variation  in  equinoctial  months  is  quite  small.  The  highest  values  (excluding  the 
less  reliable  transition  period  around  sunrise)  are  seen  in  winter  months,  around  midnight, 
at  times  of  low  solar  activity.  This  tells  us  that  the  electron  density  profiles  are  flattest 
at  these  times.  According  to  Buonsanto  et  al.  (1979)  this  occurs  to  a  greater  extent  at 
higher  latitudes  (Goose  Bay,  Canada  -  L=4).  This  all  suggests  that  the  layer  peak  has  a 
far  more  rapid  decay  than  the  TEC  after  sunset  during  winter,  at  low  activity  and  at  higher 
latitudes  and  flatter  profiles  result.  Titheridge  (1973)  suggests  that  this  increase  in  r  relates 
to  a  lowering  of  the  0+ / H+  transition  height.  The  profiles  attain  a  more  normal  shape 
(and  thickness)  before  dawn.  The  higher  slab  thickness  values  attained  during  summer 


days  could  be  interpreted  in  terms  of  greater  scale  heights  occurring  in  the  topside  region 
during  times  of  maximum  observed  TEC  (increases  in  ionization  being  more  likely  to  occur 
in  the  topside  than  the  layer  peak),  and  would  be  related  to  the  higher  electron  and  ion 
temperatures  occurring  then.  The  generally  clear  trend  to  increase  r  with  increasing  solar 
activity  (comparing  the  lower  and  upper  panels  in  Figure  2)  also  reflects  increases  in  the 
topside. 

The  next  step  is  to  determine  a  convenient  numerical  (coefficient-based) 
representation  for  the  mean  v<*.o,es  of  r  from  the  Hamilton/Waliops  Island  site.  As  the 
diurnal  variation  appeared  to  dominate,  this  was  the  first  variation  to  be  described.  Various 
numerical  techniques  were  considered  for  describing  the  diurnal  variation.  Polynomial 
fits  were  not  accurate  (up  to  a  reasonable  order)  and  are  not  intrinsically  periodic.  Cubic 
splines  are  generally  a  smooth  and  accurate  representation  of  data  but  are  more  useful 
for  interpolation  than  for  convenient  numerical  summaries,  and  can  be  unreliable  when 
individual  data  points  are  suspect.  Fourier  analysis  was  the  logical  choice  for  describing  this 
generally  smooth  and  periodic  database  (the  monthly  trimmed  means).  The  order  of  the  fit 
was  optimized  both  with  respect  to  the  residuals  of  the  fit  and  with  respect  to  the  smooth 
appearance  of  the  resulting  diurnal  curves,  and  overall,  4th  order  Fourier  analysis  yielded 
the  best  results.  Furthermore,  a  Weighted  Fourier  analysis  was  used  in  an  attempt  to 
improve  the  reliability  of  the  representation.  In  any  given  month,  weights  were  assigned  to 
each  of  the  24  hourly  mean  values  that  reflected  the  level  of  confidence  in  that  value.  These 
weights  were  inversely  proportional  to  the  quartile  range  of  the  individual  hourly  values,  so 
•hat  the  mean  of  a  widely  spread  distribution  of  values  gets  little  weight.  In  addition,  data 
obtained  during  sunrise  and  sunset  periods  were  further  downweighted  to  avoid  possibly 
erroneous  values  obtained  when  horizontal  gradients  were  large.  Weights  at  these  times 
were  further  reduced  by  a  factor  of  10.  Figure  4  contrasts  a  raw  Fourier  analysis  with  a 
Weighted  Fourier  analysis  for  the  same  set  of  data.  In  the  majority  of  cases,  there  is  little 
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difference  between  the  two.  The  only  systematic  difference  is  that  a  peak  in  mean  r  that 
is  often  seen  around  dawn  hours  does  not  appear  in  the  Weighted  Fourier  representation. 
The  justification  for  this  approach  is  that  it  would  be  difficult  to  describe  reliably  such  a 
feature  when  it  is  often  denoted  by  only  one  point  (when  hourly  data  is  used).  Further, 
to  devote  a  large  effort  to  modeling  a  peak  in  r  at  0500  LT  or  0600  LT  is  largely  pointless, 
as  the  ionosphere  is  changing  rapidly  and  the  mean  ionospheric  models  arc  all  less  reliable 
anyway.  Yet,  the  presence  of  this  dawn  peak  is  also  worthy  of  comment.  While  it  has  been 
neglected  in  our  numerical  model,  it  is  a  persistent  feature.  Indeed,  it  is  a  common  feature 
in  studies  of  slab  thickness,  having  been  observed  in  Wales  (Hajeb-Hosseinieh  and  Kersley, 
1975),  New  Zealand  (Titheridge,  1973),  Australia  (Essex,  1978)  and  India  (e.g.  Prasad  et 
al.,  1987),  as  well  as  farther  north  in  the  American  sector  (Buonsanto  et  al.,  1979).  Such 
consistency  would  certainly  indicate  a  real  feature.  Evans  (1968)  and  Titheridge  (1973) 
identified  this  feature  with  an  enhanced  loss  mechanism  near  the  layer  peak  prior  to  dawn. 
Indeed,  when  one  looks  at  individual  months  where  this  peak  is  pronounced,  the  rise  in  r 
relates  to  a  fall  in  Nmax  greater  then  the  fall  in  TEC. 

Having  described  the  diurnal  variations  by  this  Weighted  Fourier  analysis,  the 
dependencies  of  r  on  solar  and  geomagnetic  activity  were  sought.  In  particular,  the 
variation  of  r  with  respect  to  monthly  mean  sunspot  number  (a  ready  measure  of  solar 
activity)  and  monthly  mean  Ap  index  (geomagnetic  activity)  was  studied.  A  linear 
relationship  with  respect  to  sunspot  number  was  deduced.  Overall,  the  correlation 
coefficient  was  around  0.4.  While  this  indicates  significant  departures  from  a  linear 
dependence,  it  also  means  that  inclusion  of  this  simple  dependence  does  significantly  reduce 
the  overall  residuals.  There  was  no  justification  for  considering  higher  order  terms  than 
first.  It  was  noted  that  there  was  some  diurnal  modulation  in  the  degree  of  correlation 
between  r  and  sunspot  number,  with  the  greater  correlation  generally  occurring  during  the 
day  and  significantly  less  correlation  occurring  at  night,  especially  around  dawn.  Overall,  no 
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dependence  on  monthly  mean  magnetic  activity  could  be  found.  The  correlation  coefficient 
in  this  case  was  0.04.  This  was  not  a  surprising  result,  as  monthly  mean  magnetic  activity  is 
not  a  meaningful  quantity.  Day-to-day  variations  are  enhanced  during  times  of  geomagnetic 
storms  and  so  day-to-day  effects  would  be  expected,  but  it  is  not  obvious  that  these  effects 
would  carry  over  to  monthly  means.  Kersley  and  Hajeb-Hosseinieh  (1976)  did  report 
a  geomagnetic  dependence  in  the  slab  thickness  over  Wales  (this  was  also  described  in 
their  earlier  unpublished  report).  However,  this  trend  was  probably  a  chance  result.  The 
database  extended  over  less  than  one  solar  cycle  and  the  apparent  dependence  on  mean  Ap 
values  hinged  on  several  months  from  1973  and  1974  of  high  mean  Ap.  The  same  trends 
towards  higher  r  in  these  particular  months  were  in  fact  seen  in  our  Hamilton  database, 
but  no  correlation  was  seen  at  Hamilton  as  the  effect  simply  did  not  last  for  the  whole  two 
solar  cycles.  It  thus  appears  that  the  smaller  sample  size  led  to  an  apparent  correlation  with 
monthly  mean  Ap  in  the  case  of  the  Welsh  data. 

Finally,  the  annual  variations  of  the  model  were  described.  The  derived  coefficients 
of  the  diurnal  and  solar  cycle  variations  were  examined  as  a  funct  ion  of  month  (examples 
are  seen  in  Figure  5).  The  plots  indicated  that  a  ls<  order  Fourier  series  would  be  a  good 
representation  of  the  annual  variations.  There  was  no  indication  of  any  2nJ  or  higher  order 
terms.  Coefficients  not  well-defined  by  a  1st  order  analysis  were  simply  too  noisy  to  be 
described  accurately  by  any  method. 

The  coefficients  for  this  simple  mid-latitude  slab  thickness  model  are  given  in  Table 
1,  and  a  sample  Fortran  routine  to  generate  predicted  values  of  r  using  these  coefficients 
is  listed  in  the  Appendix.  The  rms  errors  of  this  model,  with  respect  to  the  original 
trimmed  means,  is  23  km,  or  less  than  10%  of  the  total  sample  average  of  277  km  (with 
corresponding  a  of  71  km).  On  a  month-by-month  basis,  the  errors  arc  normally  distributed 
(that  is,  there  is  no  bias  in  the  rms  due  to  individual  poor  months),  and  no  residual 
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dependency  on  a  monthly  parameter  could  be  detected. 

The  overall  variations  included  in  this  simple  model  are  shown  in  Figure  6  in  the  form 
of  Figure  1.  The  smoothness  and  reliability  of  the  model  can  be  gauged  by  comparing  these 
two  figures.  Figure  7  shows  the  predictions  of  the  model  of  Hajeb-Hosseinieh  and  Kersley 
(1975),  ignoring  the  probably  spurious  geomagnetic  dependence.  Apart  from  the  Hajeb- 
Hosseinieh  and  Kersley  model  predicting  slightly  higher  values  the  two  models  are  in  very 
good  agreement,  and  so  it  could  be  expected  that  at  least  over  the  40°  to  50°  range,  that  r 
is  not  a  strong  function  of  either  latitude  or  longitude. 

3.2  Day-to-day  Variations 

In  this  section  we  will  describe  both  the  results  of  a  correlation  analysis  of  the 
daily  TEC,  Nmax  and  r  values  in  our  database,  and  some  preliminary  results  from  an 
investigation  of  the  ionospheric  effects  of  geomagnetic  storms  (regarded  as  the  extreme  cases 
of  day-to-day  variability). 

We  begin  with  the  correlation  analysis.  For  this  study,  we  are  required  to  consider 
distributions  of  each  of  the  three  variables,  and  to  investigate  the  degree  to  which  a 
departure  from  the  mean  in  one  variable  is  associated  with  departures  from  the  mean  in  the 
others.  Perfectly  correlated  variables  are  distinguished  by  a  one  a  deviation  in  one  variable 
being  associated  with  a  one  a  deviation  in  the  other,  and  has  a  correlation  coefficient  of 
1.0.  A  correlation  coeficient  of  -1.0  denotes  anti-correlation,  while  0.0  would  show  that 
no  correlation  is  present  between  the  two  variables.  We  have  selected  a  sample  size  of  one 
month  for  each  hourly  value,  and  have  deduced  correlation  coefficients  for  each  month-hour 
sample  throughout  the  database.  A  histogram  of  the  correlation  coefficients  is  shown  in 
Figure  8.  As  expected,  TEC  and  foF2  (and  therefore,  Nmax)  are  quite  highly  correlated. 
The  distribution  of  values  peaks  around  0.9  and  there  are  only  a  small  number  of  samples 
with  correlation  coefficients  less  than  0.5.  Also  worth  noting,  is  the  distribution  of  values 
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for  r  and  foF2.  There  is  a  tendency  for  these  values  to  be  anti-correlated,  with  a  peak 
value  of  around  -0.7  to  -0.8.  This  anti-correlation  is  significant  in  terms  of  developing  a 
daily  ‘adjustment’  to  the  mean  slab  thickness  model.  In  a  real-time  case  a  value  of  foF'2 
may  be  available,  and  a  real-time  estimate  to  r  could  be  desired  before  predicting  TEC. 

If  the  departure  from  the  mean  foF2  value  is  used  in  combination  with  an  assumed  anti¬ 
correlation  with  r,  and  the  resulting  r  used  to  estimate  the  real-time  TEC,  this  TEC 
estimate  will  be  essentially  the  same  value  deduced  from  mean  values  of  both  foF2  and 
r.  That  is,  assuming  a  perfect  anti-correlation  tends  to  wipe  out  the  known  departure 
rather  than  to  make  use  of  it.  Similarly,  if  one  uses  ‘daily’  estimates  of  r  in  this  fashion  to 
generate  daily  values  of  TEC  from  daily  foF2,  then  the  predicted  day-to-day  variability 
in  TEC  will  be  greatly  underestimated  (that  the  anti-correlation  is  not  perfect  means 
that  the  departures  will  not  exactly  cancel).  This  does  not  mean  that  a  model  of  daily 
r  or  TEC  values  cannot  be  found,  or  that  a  description  of  the  day-to-day  variability 
cannot  be  obtained.  It  simply  means  that  a  simple  statistical  analysis  does  not  provide  the 
answers.  The  place  to  begin  a  proper  investigation  of  day-to-day  variability  is  in  the  periods 
immediately  following  geomagnetic  storms  where  the  variability  is  greatest. 

The  behavior  of  TEC  (and  therefore  Nmax,  the  two  being  well  correlated)  during 
geomagnetic  storm  periods  is  well-known.  Mendillo  (1978)  used  data  from  a  chain  of 
stations  at  American  sector  longitudes  to  study  and  compare  the  effects  at  different 
latitudes,  and  used  the  observed  variations  to  deduce  the  underlying  physical  processes. 
Mendillo  described  the  trends  with  latitude  and  season,  and  studied  some  selected  storm 
periods  in  greater  detail.  In  this  report,  some  preliminary  results  for  the  ionospheric 
response  to  geomagnetic  storms,  for  each  of  TEC,  Nmax  and  r,  from  our  two  solar  cycle 
database  are  presented.  When  studying  daily  responses,  there  are  a  couple  of  points  in 
keep  in  mind.  Firstly,  the  departures  from  mean  conditions  are  desired.  Mendillo  (1978) 
considered  the  percentage  departures  from  the  monthly  means  and  we  will  base  our  work 
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on  the  same  parameter.  Secondly,  some  care  is  required  to  investigate  the  departures 
from  normal  conditions  on  a  storm-by-storm  basis.  For  one  thing,  the  ionosphere  can  be 
sufficiently  disturbed  as  to  disrupt  the  continuity  of  data  sampling  during  storm  periods. 
Also,  actual  hourly  data  will  be  noisier  than  monthly  means  (especially  for  r)  and  will  have 
individual  spurious  values,  and  thus  some  form  of  data  screening  should  be  implemented 
prior  to  studies  of  individual  storms. 

For  now,  however,  we  simply  wish  to  present  some  results  for  two  solar  cycles 
of  Hamilton  TEC,  Nmax  and  r  storm-related  departures.  To  date,  the  data  has  not 
been  shifted  in  storm-days  in  the  manner  of  Mendillo  (1978).  Figure  9  shows  the  storm 
departures  averaged  over  all  storms  in  the  period.  The  same  trends  seen  at  Hamilton  for 
a  smaller  sample  by  Mendillo  are  seen  here.  Figures  10  and  11  show  the  effects  of  limiting 
the  storms  considered  in  these  averages  according  to  a  minimum  value  of  Ap  (a s  deduced 
from  3-hourly  Kp).  The  effects  of  increasing  Ap  are  clearly  seen,  and  thus  a  dependence  on 
Ap  for  both  storm-time  departures  and  day-to-day  variability  would  be  expected.  Figures 
12  and  13  demonstrate  the  overall  averages  (all  storms)  of  departures  during  summer  and 
winter  months.  It  was  expected  from  the  discussion  of  Mendillo  (1978),  that  the  night-time 
enhancement  in  r  would  become  more  significant  at  lower  latitudes  (in  this  case,  Hamilton) 
during  the  winter  months,  and  this  is  seen  in  the  average  winter  patterns  shown  here  for  all 
but  the  first  night. 

In  all  these  storm-departure  figures,  one  consistent  feature  is  that  r  is  enhanced  during 
times  of  geomagnetic  disturbance.  That  is,  the  electron  density  profile  remains  flatter  than 
normal,  both  in  the  TEC  enhancement  (filling  in  of  the  topside)  and  depletion  (enhanced 
loss  in  the  layer  peak)  phases.  Mendillo  (1978)  relates  the  daytime  increases  in  r  to  the 
heating  processes  that  cause  enhancements  in  the  plasma  scale  height..  Another  noticeable 
feature  in  the  response  in  t  is  that  the  variations  are  smaller  than  either  TEC  or  Nmax, 
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and  that  on  a  storm-by-storm  basis,  the  response  is  not  as  smooth  as  for  TEC  and  Nmax. 
This  reflects  the  secondary  nature  of  the  r  variations  due  to  the  generally  high  degree  of 
correlation  between  TEC  and  Nmax,  as  well  as  the  fact  that  less  data  is  available  for  r  than 
either  TEC  or  Nmax  (as  both  are  required  before  r  can  be  defined). 

It  is  important  to  recall  that  storms  are  not  all  of  the  same  character  and  responses 
vary  significantly  from  one  storm  to  the  next.  Indeed,  a  20%  average  departure  may 
correspond  to  a  f>0%  departure  for  |  of  the  cases  and  no  departure  elsewhere.  It  is  more 
useful  to  correlate  the  storms  where  the  60%  effect  occurs  than  to  model  a  20%  term  overall. 
The  next  logical  step  would  be  to  characterize  the  storms  in  terms  of  the  observed  variations 
and  to  correlate  these  with  season,  solar  activity,  maximum  Ap,  onset  time  and  so  on. 
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4.  VALIDATING  EXISTING  IONOSPHERIC  MODELS 

The  availability  of  a  mid-latitude  mean  r  model  makes  two  tasks  immediately 
possible.  The  first  is  to  convert  a  database  of  foF2  values  from  a  network  of  mid-latitude 
ionosondes  into  a  database  of  TEC  values,  to  be  used  for  model  validation.  The  second  is 
to  use  the  currently  available  profile  models  to  predict  r  at  the  location  where  it  is  known 
(Hamilton),  to  further  validate  the  models. 

The  first  of  these  tasks  was  a  project  coordinated  at  AFGL,  involving  data  converted 
from  foF2  to  TEC  at  the  Institute  for  Telecommuncation  Sciences,  Boulder,  Colorado 
using  the  simple  mean  model  of  r  developed  here  at  Boston  University,  with  personnel 
at  Computational  Physics  Inc.,  Newton,  MA,  making  comparisons  of  the  derived  TEC 
values.  This  work  is  currently  being  prepared  as  a  report  (Klobuchar  et  al.,  1989),  but 
it  is  worth  describing  a  few  key  points  here,  for  completeness.  Two  of  the  models  tested, 

IRI  (Rawer,  1981)  and  FAIM  (Anderson  et  al.,  1989),  were  consistently  poor  predictors 
of  the  deduced  TEC.  The  Bent  (Llewellyn  and  Bent,  1973),  Damon- Hartranft  (Damon 
and  Hartranft,  1970),  and  ICED  (Tascione  et  al.,  1988)  models  were  all  more  successful, 
typically  describing  the  TEC  to  around  one  a  of  the  apparent  day-to-day  variation.  It 
was  noted  that  a  significant  portion  of  the  errors  in  the  deduced  TEC  could  be  related  to 
the  predicted  foF2  in  each  model  not  always  being  in  good  agreement  with  the  observed 
monthly  median  foF2. 

Here  at  Boston  University,  available  empirical  electron  density  profile  models  have 
been  used  to  predict  the  mean  values  of  r  at  Hamilton,  MA.  Some  comparisons  showed  that 
none  of  the  models  tested  had  r  values  that  were  strong  functions  of  latitude  or  longitude 
in  this  mid-latitude  range,  and  therefore  these  comparisons  are  equally  valid  at  other  mid- 
latitude  locations.  Estimates  of  r  have  been  generated  from  each  of  the  IRI,  FAIM,  Bent 
and  Damon-Hartranft  models,  and  3-D  plots  of  these  are  shown  in  Figures  14-17.  As  was 
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the  case  for  the  TEC  comparisons,  the  IRI  and  FAIM  models  are  an  inadequate  description 
of  r.  In  the  IRI,  there  is  little  diurnal  or  seasonal  variation,  and  an  exaggerated  variation 
with  solar  activity.  FAIM  showed  some  diurnal  variation,  but  little  seasonal  and  solar  cycle 
variation.  The  Bent  model  r  values  were  of  consistently  different  character  to  the  observed 
values,  and  so  it  was  a  little  surprising  that  the  derived  TEC  values  should  have  been  quite 
comparable.  The  Damon-Hartranft  model  showed  the  intrinsic  variation  of  r  that  was  most 
like  the  Hamilton  observed  values,  but  contained  somewhat  larger  diurnal  modulation. 

As  both  the  Bent  and  Damon-Hartranft  profile  models  were  close  overall  in  TEC,  and 
of  simple  enough  formulation  to  be  changed,  some  attempts  were  made  to  adjust  each  of 
these  models  to  match  the  observed  Hamilton  r. 

The  Bent  model  profile  consists  of  five  regions,  and  is  described  by  a  bottomside  bi¬ 
parabola,  topside  parabola,  and  three  topside  exponential  layers.  The  parabolic  layers  both 
yielded  slab  thickness  estimates  that  were  proportional  to  the  layer  semi-thickness,  ym,  and 
these  were  largely  responsible  for  the  daytime  variations  in  r.  The  first  adjustment  made 
was  therefore  to  replace  the  tabulated  values  of  ym  with  values  that  were  proportional  to 
observed  values  at  Hamilton.  This  improved  the  agreement  in  the  afternoon  and  evening. 
Nighttime  contributions  were  dominated  by  the  topside  exponentials,  and  large  values 
between  midnight  and  dawn  could  be  remedied  by  reducing  the  scale  height  of  these 
exponentials  in  certain  cases.  The  adjusted  Bent  model  r  values  are  shown  in  Figure  IS. 
While  agreement  is  still  not  ideal  (both  because  the  variations  are  still  quite  different  and 
because  of  the  circularity  of  implicitly  using  observed  values  to  reproduce  observed  values) 
there  is  considerable  improvement  over  the  original  Bent  model. 

The  Damon-Hartranft  model  consists  of  only  four  layers;  E,  FI,  bottomside  F2  and 
topside  F2.  The  first  three  are  Chapman-type  layers,  while  the  topside  F2  is  described  by 
an  exponential  layer  of  increasing  scale  height.  In  terms  of  r,  the  contribution  from  the 
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topside  dominated.  However,  there  was  an  unexpectedly  large  contribution  from  the  Fl 
layer  (over  20%  in  some  cases)  that  peaked  during  the  afternoon  when  the  predicted  value 
of  r  was  so  much  higher  than  that  observed.  Without  this  Fl  layer,  the  Damon-Hartranft 
model  r  values  are  as  shown  in  Figure  19.  This  is  the  best  agreement  of  any  ’model’  with 
the  Hamilton  slab  thickness,  and  would  therefore  be  the  first  to  be  recommended  for  future 
use  based  on  these  comparisons. 

5.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  behavior  of  the  equivalent  slab  thicknes  has  been  investigated  at  a  mid-latitude 
site,  and  had  been  described  by  a  simple  numerical  model.  The  usefulness  of  this  model  at 
other  locations  depends  on  the  slowly-varying  nature  of  r,  as  deduced  from  comparisons 
made  either  side  of  the  Atlantic  Ocean.  Given  the  desireability  of  having  a  global  model  of 
r,  both  physically  and  operationally,  it  would  be  useful  to  further  investigate  the  variations 
of  r  at  other  locations.  It  would  also  be  desireable  to  obtain  TEC  measurements  from  other 
independent  sources,  as  both  the  Hamilton  and  Wales  TEC  data  was  validated  in  the  same 
way;  relying  on  comparisons  of  derived  r  values  to  correctly  identify  the  nn  ambiguity. 
Removing  any  possible  circularity  of  this  procedure  would  be  a  major  step  in  correctly 
defining  the  real  variations  of  slab  thickness. 

Deriving  the  real  day-to-day  variations  also  remains  a  desireable  goal.  Characterizing 
the  responses  following  geomagnetic  storms  will  be  a  key  step  in  pinpointing  the  physical 
processes  that  are  occurring,  as  well  as  making  possible  some  functional  model  of  day-to-day 
variations.  It  is  through  a  method  such  as  this  that  the  most  reliable  real-time  estimates  of 
TEC  (or  range  error)  will  be  achieved. 
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Figured.  Diurnal  variations  of  moan  slab  1  hickness  h.r  each  nmnlliof  I?t7l.  with  the 
months  as  labeled.  Six  line  styles  have  been  used,  each  style  denoting  two  months  spaced  six 
months  apart  (e.g.  January  and  July). 
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TABLE  1  . 


Coefficients  of  the  numerical  model 
of  slab  thickness  over  Hamilton,  MA 


967.669128 

-55.940933 
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APPENDIX. 


An  example  of  a  Fortran  routine  that  evaluates  the  monthly  mean 
slab  thickness  using  the  listed  coefficients.  Takes  month,  Local  Time  and 
monthly  sunspot  number  as  input 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccca.cccccccccccccc 

program  maketau 
c 

c  Use  the  file  of  coefficients  to  predict  slab  thickness 
c  for  a  specified  month,  hour  and  sunspot  number, 

c 

integer  orderd,  order a 

real  aazero(9),  aan(4,9),  abn(4,9) 
real  bazero(9),  ban(4,9),  bbn(4,9) 

rea  I  a (9)  ,  b (9) 

real  coeffs(9) 

real  azero,  an(4),  bn(4) 

integer  month,  hour 
real  ssn 

real  sines(9,24),  coses(9,24) 
real  as i nes (9 , 12) ,  acoses(9,12) 


c 

c  Input  the  month,  ssn,  and  hour  (LT) 
c 

write(6,*)’  Enter  the  desired  month,  ssn  and  Local  Time’ 
read (5,*)  month,  ssn,  hour 

c 

c  set  up  the  convenient  array  of  cos  and  sin  values  for  the  Fourier  fits 
c 

p i  =  acos(-l .0) 
scale  =  pi /12.0 
asca I e  =  p i /6 . 0 

do  101  i  =  1,  9 
c  diurnal  (hours  LT  :  00-23) 
do  102  j  =  1,  24 

angle  =  float(i  *  (j-1))  *  scale 
coses(i,j)  =  cos(angle) 

102  sines(i,j)  =  sin(angle) 
c  and  annual  (months  1-12) 

do  103  j  =  1,  12 

angle  =  float(i  *  j)  *  ascale 

acoses(i,j)  =  cos(angle) 

103  asines(i,j)  =  sin(angle) 

101  continue 

c 

c  Read  in  the  coefficients 
c 

open (un i t=l , name= ’coef  f  s . dat ’ , status= ’o I d ’ , f  orm= ' f  or mat ted ’) 


ID 


c  define  the  orders  of  the  Fourier  fits  (Diurnal  and  Annual) 
orderd  =  4 
ordera  =  1 

c  read  in  those  annual  coeffs 

do  200  k  =  1,  2*orderd+l 

read(l,*)  j,  aazero(k) ,  (aan(j , k) , j=l .ordera) ,  (abn (j , k) , j =1 , ordera) 
read(l,*)  j,  bazero(k),  (ban (j , k ), j =1 , ordera) ,  (bbn(j , k) , j=l .ordera) 
200  continue 

c lose(l) 

c 

c  First,  use  the  annual  Fourier  stuff  to  generate  the  coefficients  with 
c  respect  to  ssn  for  the  desired  month 
c 

do  1000  k  =  1,  2*orderd+l 
tmpa  =  aazero(k)/2 .0 
tmpb  =  bazero(k)/2.0 
do  1200  n  =  1,  ordera 

tmpa  =  tmpa  +  aan(n.k)  *  acoses (n .month)  +  abn(n.k)  *  as i nes (n , month) 
tmpb  =  tmpb  +  ban(n.k)  *  acoses (n , month)  +  bbn(n.k)  *  as i nes (n , month) 
1200  continue 

a(k)  =  tmpa 
b(k)  =  tmpb 
1100  continue 

1000  continue 

c 

c  Then,  determine  the  diurnal  Fourier  coeffs  for  that  ssn 

c 

do  4000  k  =  1,  2*orderd+l 
4000  coeffs (k)  =  a(k)  +  b(k)  *  ssn 

c  or,  in  the  more  usual  Fourier  parlance 

azero  =  coeffs(l) 
do  4010  k  =  1,  orderd 
an(k)  =  coeffs(k+l) 

4010  bn(k)  =  coef f s (k+orderd+1) 

c 

c  Finally,  regenerate  the  hourly  value  using  those  diurnal  coeffs 
c 

tau  =  azero/2.0 
do  2410  n  =  1,  orderd 

2410  tau  =  tau  +  an(n)  *  coses(n,hour  +  l)  +  bn(n)  *  s i nes (n , hour* 1) 

write(6,*)’  Thus  Tau  =  tau 

stop 

end 
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